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copending U.S. application No. 10/136,813, filed April 30, 2002, entitled METHOD 
AND APPARATUS FOR IMPROVING IMAGE CLARITY AND SENSITIVITY IN 
OPTICAL COHERENCE TOMOGRAPHY IISENG" DYNAMIC FEEDBACK TO 
CONTROL FOCAL PROPERTIES AND COHERENCE GATING, both commonly 
5 assigned to the assignee of the present application, the disclosures of which are 
incorporated by reference in its entirety herem. 

FIELD OF THE INVENTION 

The present invention relates to apparatus, method, logic arrangement and storage 
medium for dramatically increasing the sensitivity in tiie detection of ojptical 
10 coherence tomography and low coherence interferometry signals by detecting a 
parallel set of spectral bands, each band bemg a unique combination of optical 
firequencies. 

BACKGROUND OF THE ART 

Two methods cim'entiy exist to implement depth ranging in turbid media. The first 
15 method is known as Low Coherence Interferometry ("LCI"). This method uses a 
scanning system to vary the reference arm length and acquire the interfiarence signal at 
a detector and demodulating the fiinge patt^ to obtain the coherence envelope of the 
source cross correlation function. Optical coherence tomography ("OCT") is a means 
for obtaining a two-dimensional image using LCL OCT is described by Swanson et al. 
20 m U.S. Pat. No. 5,321 ,501 . Multiple variations on OCT have been patented, but many 
suffer firom less than optimal signal to noise ratio ("SNR"), resulting in non-optimal 
resolution, low imaging fi-ame rates, and poor depth of penetration. Power usage is a 
&ctor in such imaging techniques. For example in ophthalmic uses, only a certain 
number of milliwatts of power is tolerable before thermal damage can occur. Thus, 
25 boosting power is not feasible to inorease SNR in such environments. It would be 
desirable to have a method of raising the SNR without appreciably increasing power 
requirements. 

A second method for depth ranging in turbid media is known in the literature as 
spectral radar. In spectral radar the real part of the cross spectral density of sample and 
30 reference arm light is measured with a spectrometer. Depth profile information is 
encoded on the cross-spectral density niodulation. Prior designs for spectral radar is 
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primarily found in the literature. 

The use of spectral radar concepts to increase the signal to noise ratio of LCI and OCT 
have been desaibed earlier. However, in this description, only the real part of the 
• complex spectral density is measured and the method uses a large number of detector 
5 elements (about 2,000) to reach scan ranges on the order of a millimeter. It would be 
desirable to have a method that would allow for an arbitrary number of detector 
elements. Secondly, the previously described method uses a single charge coupled 
device ("CCD") to acquire the data. Since the charge storage capacity is limited, it 
requires a reduction of the reference arm power to approximately the same level as the 
10 sample arm power, giving rise to auto correlation noise on the sample arm light. In 
addition, since no carrier is generated, the 1/f noise will dominate the noise in this 
system. Thirdly, even with the short integration times of state of the art CCD 
technology, phase instabilities in the interferometer reduce fringe visibility of the 
cross spectral density modulation. 

15 SUMMARY OF THE INVENTION 

The present invention can increase the SNR of LCI and OCT by splitting the LCI 
broad bandwidth source into a number **N" of spectral bands. In one exemplary 
embodiment, the N spectral bands are individually detected and processed to provide 
an increase in the SNR by a factor of N. This increase in SNR enables LCI or OCT 

20 imaging by a factor of N times faster, or altematively allows imaging at tihe same 
speed with a source that has N times lower power. As a result, the present invention 
overcomes two of the most important shortcomings of conventional LCI and OCT, 
namely, source availability and scan speed. The factor N may reach more than 1,000, 
and allows construction of OCT and LCI systems that can be more than three orders 

25 of magnitude improved from OCT and LCI technology currently in practice. 

The present invention improves current data acquisition speeds and availability of 
sources for OCT. Shot noise is due to flie statistical fluctuations of the current that are 
due to the quantized or discrete electric charges. The reduction of shot noise allows 
for much lower source powers or much higher acquisition rates. Limitations in current 
30 data acquisition rates (approximately 4 frames/sec) are imposed by available source 
power and availability of fast mechanisms for scanning delay. An increase in the 
sensitivity of the detection by a factor of 8 would allow real time imaging at a speed 
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of about 30 frames per second. An inorease of the sensitivity by a factor of about 
1,000-2,000 would allow for the use of sources with much lower powers and higher 
spectral bandwidths which are readily available, cheaper to produce, and can generate 
higher resolution LCI or OCT scans. 

5 For ophtihalmic applications of OCT, the efficient detection preferably allows for a 
significant increase of acquisition speed. The limitation in ophthalmic applications is 
the power that is allowed to enter tiie eye according to the ANSI standards 
(approximately 700 microwatts at 830 nm). Current data acquisition speed in 
ophthahnic applications is approximately 100-500 A-lines per second. The power 
10 efficient detection would allow for A-Une acquisition rates on the order of about 
100,000 A-lines per second, or video rate imaging at about 3,000 A-lines per image. 

The gain in SNR is achieved because the shot noise has a white noise spectrum. An 
intensity present at the detector at frequency co (or wavelength X) contributes only to 
the signal at frequency cl>^ but the shot noise is generated at aU frequencies. By 
15 narrowing the optical band width per detector, the shot noise contribution at each 
frequency can be reduced, while the signal component remains the same. 

In summary, the present invention improves a performance of LCI and OCT, and as a 
result, can be used in developing LCI and OCT diagnostic technologies for medical 
and non-medical applications. 

20 Other features and advantages of the present invention will become apparent upon 
reading flie following detailed description of embodiments of the invention, when 
taken in conjunction with the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is illustrated in the drawings in which like reference characters 
25 designate the same or similar parts throughout the figures of which: 

Fig. 1 is a schematic view of a conventional system. 

Fig^ is a schematic view of a preferred embodiment of the parallel detection scheme 
forLCL 

Fig^ is a schematic view of a system with one detector array according to one 
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embodiment of the present invention. 
F ig. 4 i s a detail of a probe. 

Fig. 5 is a schematic view of separating unit in combination with two integrating CCD 
arrays for detection of the dual-balanced wavelength demultiplexed signal. 

5 Fig. 6 is a schematic view of a preferred embodiment of a standalone system 

Fig/T^is a schematic view showing spectral separating into 2 bands. 

Fig. 8 is a schematic of spectral separating into 4 bands. The spectral resolution 
preferably used for each detector is twice as coarse as in the case of multiplexing into 
2 bands. 

10 Fig. 9 is a schematic view of using beam recombination to provide one dimension of 
interference information along one dimension of a two-dimensional detector array, 
while performing wavelength separating along the other dimension of the two 
dimensional array. 

Fig. 10 is a schematic view of a phase tracking system according to one embodiment 
15 of the present invention. 

Fig. 11 is a flowchart depicting the reconstruction of LCI or OCT signal from 
wavelength bands. 

Ijg ^l2 i s a schematic view of a spectral domain OCT interferometer design with a 
source combining the spectra of several superluminescent sources. 

20 Fi g. 13 is a schematic view of a system with a four detector array. 

Fi g^ 14 i s a graph of a typical interference patter as a ftmction of path length 
difference between sample arm and reference arm. 

Fig» 15 is an embodiment of a phase tracker system with an extended phase lock 
range. 

25 Figs. 15A-C are flow diagrams of a method. 
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FigJ4 is a graph of ftequency versus OCT power spectrum. 

Fig. 17 is a graph of firequency versus amplitude spectrum subtracted from the shot 
noise (experimental data) for the N-1 (dotted line) and N~l/3 (solid line) cases. 

Fig. 18 is a graph of power density for the full spectrum as a function of frequency. 

5 Fig. 19 is a graph after subtraction of the shot noise levels. 

Fig. 20 is a graph after processing the signals. 

Fig. 21 is a graph of the coherence envelope for the coherently summed channels. 
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Overview 

1 0 Certain exemplary onbodiments of the present invention include a hybrid method that 
implements aspects of LCI and OCT where the ref^ence aim is scanned, and spectral 
radar, which does not require reference arm scanning. 

In one embodiment, the signal in the detection arm of an OCT system is split into 
more than one spectral band before detection. Each spectral band is detected by a 

15 separate photodetector and amplified. For each spectral band, the signal can be band 
pass filtered around the signal band by analog electronics and digitized, or, 
alternatively, the signal may be digitized and band pass filtered in software. As a 
consequence, the shot noise contribution fo the signal can be reduced by a factor equal 
to the nxraiber of spectral bands, while output of the signal remains the same. The 

20 reduction of the shot noise increases the dynamic range and sensitivity of the system. 

In another exemplary embodiment of the present invention, an apparatus is provided 
for spectral radar that does not require reference arm scanning. For many detectors^ no 
ranging or reference arm scanniag is needed, and the method may be similar to the 
method which can be employed for a spectral radar except tiiat phase information of 
25 the cross spectral density is preferably preserved. 

In other exemplary embodiments, the present invention describes an arrangement for 
spectral radar that eliminate phase instability in the interferometer, obtaining the 
complex spectral density and eliminating auto correlation noise on the sample aim 
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light, relative intensity noise, and 1/f noise. 
THEORY 

TIME DOMAIN VERSUS SPECTRAL DOMAIN OCT 

Nearly all conventional OCT systems are based on Time Domain scanning. In such 
5 conventional systems, the length of the reference arm in a Michelson interferometer is 
rapidly scanned over a distance corresponding to the imaging depth range. An 
alternative procedure to scanning the reference arm, is one that measures the cross- 
spectral density at the detection arm of the Michelson interferometer using a 
spectrometer. In Spectral Domain OCT, no mechanical (e.g., motionless) scanning of 
10 fhe reference arm is required, while an apparatus for generating a phase shift can be 
used. Only recently was it recognized that a significant signal to noise gain can be 
achieved by direct measurement of the cross-spectral density. 

Fig. 1 shows a schematic of a conventional Time Domain OCT system. On scanning 
the reference arm path length, interference fi±ages are formed corresponding to 
15 positions that match the distance to the three structures drawn in the sample volimie. 
A single detector is used to detect the interference fringes. By envelope detection of 
the fringe patterns, an image is constructed that maps tissue reflectivity to a given 
locations 

Certain exemplary embodiments of the present invention provide a detection principle 
20 based on Spectral Radar concepts (further referred to as Spectral Domain OCT) or a 
hybrid method between Spectral Domain and Time Domain OCT that can be more 
sensitive than current state of the art Time Domain OCT, allowing a substantial 
increase in the acquisition speed to resolution ratio. 

Principle of shot noise reduction in Spectral Domain OCT 

25 The best signal to noise performance of Time Domain OCT systems is obtained when 
the noise is shot noise limited. Shot noise can be reduced significantly by replacing 
the single element detector with a multi-element array detector. When the detection 
arm light is spectrally dispersed on the array detector, each element of the array 
detects a small wavelength fraction of the spectral width of the source. The shot noise 

30 is preferably reduced by a factor equal to the number of elements of the array. The 
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principle of the signal to noise improvement is based on the white noise characteristic 
of shot noise and the observation that only electromagnetic waves of the same 
wavelength produce interference firinges. 



5 proportional to the current (or equivalently the optical power times the quantum 
efficiency) generated in the detector. For a monochromatic beam of wavelength Xj 
entering the interferometer, the fringe frequency or carrier / at the detector is 
determined by the velocity v of the mirror, = 2v/Xj. The shot noise is proportional to 
the power (or spectral density S((o)) at wavelengfli Xj. A second wavelength is 

10 preferably coupled into the interferometer. A second fringe frequency or carrier at 
frequency^ = 2v/l^ is simultaneously present The shot noise at this second frequency 
is preferably the sum of the shot noise generated by the optical power at wavelength A; 
and X^. Also, at frequency the shot noise is the sum of the shot noise generated by 
the optical power at wavelength Xj and X2. Thus, at both frequencies a cross-shot noise 

15 term is generated by the simultaneous presence of both wavelengths at the detector. 
By ^ectrally dispersing each wavelength to a separate detector, the cross shot noise 
term can be eliminated. In this way. Spectral Domain OCT offers a significant 
improvement of signal to noise ratio over Time Domain OCT systems. 

Signal to noise analvsifg nf Time Domain versus Spectral Domain OCT. 

20 Signal 

Analysis of the Signal to Noise Ratio (SNR) in Time Domain OCT has been 
described in related publications. The interference fringe peak amplitude in time 
domain OCT is given by 



25 with Pref^Psample reference and sample arm power in Watts, respectively. In 
terms of electrical power at the detector, the signal in units [-4^] is defined as 



with Tj the quantum efficiency, ethe charge quantum and Ey^hclXihst photon 



The shot noise power density N^j^^,{f) (in units [W/Hz],[a^ /Hz] or //fz]) is 




(1) 



(2) 
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energy. The reference and sample arm powers are given by the respective reflected 
spectral densities, 

Pref^a^nple = J ^rtf^a^nplei^)^^ • (3) 

Assxraiing that the reference and sample spectral densities are equal to the source 
5 spectral density 15(0), where the sample arm spectral density is attenuated by a large 
factor, i.e., S„y (co) = S{a}\ 5,^^,^ {m) = a S{a)) with a « 1 , and inserting the above 
expression of reference and sample arm into the original definition of the signal gives, 

S = 7j^e^a]j S{a))dw] I El. (4) 

Thermal, shot noise and relative int^isitv noise contributions 

10 Three contributions to the total noise of OCT signals are: thermal noise, shot noise 
and relative intensity noise. Themial noise is generated by the feedback resistor, shot 
noise is related to the finite nature of the charge quantum resulting in statistical 
fluctuations on the current, and relative intensity noise is related to the temporal 
fluctuations due to chaotic character of classical light sources. These three 

15 contributions to the noise density in units [A^ /Hz] are given by, 

^^(/) = f^.^.2(^)V^. (5) 

k is Boltzmann's constant, T the temperature in Kelvin, Rji, the value of the 
feedback resistor, and r^^ij the coherence time of the source. Coherence time is 
related to the full spectral width at half maximum AA of a Gaussian source by the 
20 following relation, r^,, = V2ln27^;iJ/(c A>i). Shot noise limited detection is 
achieved when the second term in Eq. (5) dominates the other noise contributions. 

Signal to Noise Ratio (SNK) 

The signal to noise ratio (SNR) is given by 
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SNR^ ^ , (6) 

with BW the signal bandwidth, and parameters S and N„otseif) as described above. 
Space and frequency domain description of the OCT signal 

The OCT signal is most easily described in the space domain. For a single object in 
5 the sample arm, the interference term of the OCT. signal is proportional to the real part 
of the Fourier transform of the source spectrum S{a)) , 

/(Az) oc ReJ exp(iMz)5(A:) dk , (7) 

with Az the path length difference between sample and reference arm and k the wave 
vector. As a function of time, the OCT signal is given by, 

1 0 l{t) oc ReJ expiliotv I c)S{(xi) do) , (8) 

with v the reference arm mirror velocity. The frequency spectrum of the signal is 
given by a Fourier transform of the signal in the time domain, resulting in a complex 
function. The absolute value of this function is equal to the spectral density, 

|/(/)| = |J/(0^''"^'rf^ - SiTvfcIv), (9) 

15 which shows that the signal bandwidth is directiy proportional to the source spectral 
width and scales linearly with the reference arm mirror velocity, i.e., imaging speed. 
Eq. (9) also directiy relates the absolute value of the frequency spectrum, |/(/)|, to 
the signal S (Eq. (4)). 

Eq (9) also demonstrates that each angular frequency of the light source or 
20 equivalentiy each wavelength of the source is represented at its own frequency in the 
measured interferometric signal. The depth profile information I(t) can be obtained 
from the complex cross spectral density /(/) by a Fourier transform. 

The complex cross spectral density can also be obtained by splitting the signal I(t) in 
several spectral bands using a dispersive or interferometric element. At each detector, 

10 



wo 03/062802 



PCT/US03/02349 



only part of the complex cross spectral density is determined. Combining the cross 
spectral densities of each detector, the full spectral density of the signal is retrieved. 

Thus, the same information can be obtained by separating spectral components to 
individual detectors. Combining the signal of all detectors in software or hardware 
5 would result in the same signal as obtained with a single detector. 

Signal to noise gain with Spectral domain OCT 

In the detection arm, the spectrum can be split into two equal halves, where two 
detectors each detect one half of the spectrum. According to Eq (9), the jfrequency 
spectra at detectors 1 and 2 are given by |/i(/^ = 5(;^/v)for /</o, /,(/)=0 for 

10 />/oand lAf) = 0 for f < W\l2{f} = S{nfc/v) for />/o, respectively. The 
frequency spectrum as would be acquired by a single detector in time domain OCT is 
given by the sum of I,{f) mdl^if); /(/)=/,(/)+ /^i/)- Thus, the signal S after 
combining the spectra is equal, However /j(/) = 0 for / > /o and 12(f)- 0 for 
/ < /o , the bandwidth BW p&r detector can be reduced by a factor of 2, 

15 The noise is determined by the sum of the shot noise contributions at detectors one 
and two. From Eqs. (5) and (6), the shot noise per detector is proportional to the 
reference arm power at the detector times the bandwidth for the detector. Since the 
spectrum was split in equal halves, the reference power at detectors 1 and 2 is, 
respectively, 

20 i'^=0.5P^. P^-'O.SP^. (10) 

The sum of the shot noise contribution for the two detectors is, 

NnL <^ Pr^ X + i>4^ X 0.5BW = O.SP^^BW , (1 1) 

which may compared with the shot noise of a single detector in time domain OCT, 

NZe^Pr^BW. (12) 

25 Thus, by spectrally dispersing the detection arm light over two separate detectors, the 
signal remains the same, while the noise is reduced by a factor of 2, resulting in a net 
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SNR gain by a factor of 2. 

Extending the above analysis, it can be demonstrated that the shot noise contribution 
is reduced by a &ctor equal to the number of detectors. The sum of shot noises for N 
detector elements, where each detector element receives one N*** of the total reference 
5 power, is given by, 

noise E, N 

The signal is the same as in Time Domain OCT, and the SNR ratio for Spectral 
Domain OCT is given by. 



^noise ^E^BW 

10 Thus Spectral Domain OCT enables a SNR improvement over Time Domain OCT of 
a hundred to a thousand fold, depending on the number of detector elements N. Using 
a charge coupled array or an integrating device as a detector, such as, but not limited 
to, a line scan camera, the ratio N/BW is replaced by the integration time of the 
array, which results in, 

15 (15) 
Advantages 

The exemplary embodiment of the preset invention reduce shot noise and other 
forms of noise which allows for much lower source powers, or much higher 
acquisition rates than current systems. The increased detection sensitivity allows for 

20 real time imaging. Such imaging speed can help practitioners where motion artifacts 
are a continuing problem, such as in gastrointestinal, ophthalmic and art^al imaging 
environments. By increasing the frame rate while maintaining or improving the signal 
to noise ratio such artifacts can be minimized. The present invention also enable one 
to screen large areas of tissues with OCT and allows clinical viable screening 

25 protocols using this method. 
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Fig. 2 shows a top level system 100 configiiration according to an exemplary 
embodiment of the present invention, which comprises an interferometer 102 with a 
source arm 104, a sample arm 106, a reference arm 108, and a detection arm 110 with 
a spectral separating vmit 112, multiple detectors 114, amplifiers 116, optional analog 
5 processing electronics 118 (not shown, but known to those skilled in the art), and A/D 
converters 120 (not shown, but known to those skilled in the art) for conversion of 
signals. A processing and display unit 122 has optionally digital band pass filtering 
C*BPF") units 124, Digital Fast Fourier Transforais ("FFTs") 126 (not shown), 
coherent combination of signals, and data processing and display algorithms. The 

10 detector array 114 may be IxN for simple intensity ranging and imaging and/or 
Doppler sensitive detection, 2xN for dual balanced detection, 2xN for simple intensity 
ranging and/or polarization and/or Doppler sensitive detection, or 4xN for combined 
dual balanced and polarization and/or Doppler sensitive detection. Alternatively, an 
MxN array may be used for arbitrary number **M" of detectors 1 14 to allow detection 

IS of transverse spatial information on the sample 130. 

Fig. 3 shows a schematic of one exemplary embodiment of a Spectral Domain OCT 
system 200, which includes a light source 202, splitter 204, reference arm 206, sample 
arm 208, tissue sample 130, optical element 210, grating 212, lens 214, detector 216 
array, and processor 218. The detection arm light is dispersed by the grating 212 and 
20 the spectrum imaged onto a detector array 216. By stepping the reference amf'206 
lengtih over a distance X/8, the cross spectral density of reference arm 206 and sample 
aim 208 ligfht can be determined. A Fourier transform of the* cross spectral density 
generates the depth profile information. 

Sources 

25 The source arm 203 contains at least light source 202 that is used to illimiinate the 
interferometer with low-coherence light. The source temporal coherence length is 
preferably shorter than a few microns (a preferred range is about 0.5 jam - 30 jim). 
Examples of sources include, but are not limited to, semiconductor optical amplifier, 
superluminescent diodes, light-emitting diodes, solid-state femtosecond sources, 

30 amplified spontaneous emission, continuum sources, thermal sources, combinations 
thereof and the like. Other appropriate sources known to those skilled in the art may 
be used. While ligjit is referred to herein as the source, it is intended that other 
electromagnetic radiation ranges may be suitable for use, depending on the 
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circumstances* 
Interferoineter 

The sample ami 208 collects light reflected from the tissue sample 130 and is 
combined with the light from the reference arm 206 to form interference fringes. The 
5 reference arm 206 returns light back to be combined with the source arm 203. The 
reference arm can also be transmissive with no reflection. This action of beam 
splitting/recombining may be performed using a beam splitter 204 (Michelson), or 
circulator(s) (Mach-Zehnder) or otiier means known to ttiose skilled in the art for 
separating a beam into multiple paths and recombining these multiple beams in a 
10 manner that interference between tiie beams may be detected. The splitting may be 
accomplished in free space or by using a splitter 204 having passive fiber optic or 
waveguide components. 

Sample arm 

For LCI applications, the sample arm may be terminated by an optical probe 

15 comprising a cleaved (angled, flat, or polished) optical fiber or free space beam. A 
lens (such as, but not limited to, aspherical, gradient index, spherical, diflfractive, ball, 
drum or the like) may be used to focus the beam on or within the sample. Beam 
directing elements (such as, but not limited to, mirror, prism, diffractive optical 
element or the like) may also be contained within the probe to direct the focused beam 

20 to a desired position on the sample. For OCT applications, the position of the beam 
may be changed on the sample as a fiinction of time, allowing reconstruction of a two- 
dimensional image. Altering the position of tiie focused beam on the sample may be 
accomplished by a scanning mirror (such as, but not limited to, a galvanometer, 
piezoelectric actuator or the like), electrooptic actuator, or moving the optical fiber 

25 (for example, rotating the optical fiber, or linearly translating the optical fiber). The 
sample arm probe may be a fiber optic probe that has an internally moving element 
where the motion is initiated at a proximal end of tiie probe and the motion is 
conveyed by a motion transducing arrangement (such as, but not limited to, wire, 
guidewire, speedometer cable, spring, optical fiber and the like) to the distal end. The 

30 fiber optic probe may be enclosed in a stationary sheath which is optically transparent 
where the light exits the probe at the distal end. Fig. 4 shows a detail view having an 
iimer cable 260 (which may rotate or linearly translate along the axis of the probe), an 
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outer transparent or semi-transparent sheath 262, distal optics 264, and remitted light 
266 (which may be at any angle with respect to axis of catheter). 

Reference arm delay 

A mechanism 270 in the reference arm 206 allows for scanning the group delay of the 
5 reference arm 206, This group delay can be produced by any of a number of 
techniques known to those skilled in the art, such as, but not limited to, stretching an 
optical fiber, free space translational scanning using a piezoelectric transducer, or via 
a grating based pulse shaping optical delay line. Preferably, the delay is introduced by 
a non-mechanical or motionless arrangement By ^'non-mechanical" it is meant that no 

10 mechanically moving parts are utilized. The absence of mechanically moving parts is 
believed to reduce the known deficiencies of using mechanical devices to introduce 
delay. As opposed to traditional LCI or OCT systems described in the literature, the 
reference arm 206 in the present invention does not necessarily need to scan over the 
full ranging depth in the sample, and preferably scans over at least a fraction of the 

15 ranging depth equal to one over the number of detectors (1/N)- This scanning feature 
is fundamentally different from known delay scanning schemes used in convmtional 
known LCI and OCT systems. The reference arm 206 optionally has a phase 
modulator mechanism (described more fiilly herein), such as, but not limited to, an 
acoustooptic modulator, electrooptic phase modulator or the like, for generating a 

20 carrier frequency. In order to reduce the scan range of the reference arm 206, the 
spectrum is preferably spUt into a plurality of spectral bands according to a method 
that will be explained below. 

Detection 

Referring to Fig. 2, in the detection arm 110 spectral separating unit separates the 
25 spectral components and the signal is forwarded to separate detectors 114. The 
detectors 114 may preferably consist of photodiodes (such as, but not limited to, 
silicon, InGaAs, extended InGaAs, and the like). Alternatively, a one or two 
dimensional array of detectors 1 14 (such as, but not limited to, photodiode array, 
CCD, CMOS array, active CMOS array, CMOS "smart pixel" arrays, combinations 
30 thereof and the like) may be employed for detection. Two detectors 114 for each 
spectral band may be used for polarization sensitive detection following separation of 
the recombined light into orthogonal polarization eigenstates. Detector 114 arrays 
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may be IxN for simple intensity ranging and imaging and/or Doppler sensitive 
detection, 2xN for dual balanced detection, 2xN for intensity ranging and imaging 
and/or polarization sensitive and/or Doppler sensitive detection, or 4xN for combined 
dual balanced and intensity ranging and/or Doppler sensitive and/or polarization 
5 sensitive detection. Alternatively, an MxN array may be used for arbitrary M to allow 
detection of transverse spatial information on the sample 40. 

Detector signals can be amplified by Trans Impedance Amplifiers ("TIAs") 116, band 
pass filters 124 (digitally or using analog circuitry) and digitized by AID converters 
and stored in a computer 122 for further processing. Each detector 1 14 is preferably 
10 configured to be shot noise limited. Shot noise limited detection is preferably 
achieved by adjusting the intensity of light retumed firom the reference arm 108 so 
' that the shot noise dominates over the themial noise of the resistor in the TIA 116 and 
is higher flian the relative intensity noise (**RIN"). Each detector 114 is balanced for 
such dual noise reduction^ 

15 In one embodiment of the present invention, the number of detectors 1 14, N can be in 
the range of 2-10,000 or more. A preferred range of N is about 8-10,000 detectors. In 
one preferred embodiment, eigjit detectors 1 14 (or a number in that area) can provide 
real time, or close to real time, imaging. 

Alternatively, another way for detection includes an integrating one-dimensional or 
20 two-dimensional detector 1 14 array which is capable of obtaining images at a rate 
preferably greater than 1/f noise (f = firequency) (see Fig. 5). Optionally, the BPF can 
be implemented discretely following digitization. An additional modification includes 
using an optional second detector 115 array for balanced detection which allows 
increased reference arm power and acquisition speed due to reduction of RIN and 1/f 
25 noise. In a preferred embodiment, a phase tracking apparatus and/or algorithm is used 
in the reference arm 108 to reduce signal attenuation due to fiinge instability. 

This system could be implemented using a single detector 114 with dual-balanced 
detection enabled by either interleaving dual balanced rows of the array detector or by 
placing two similar array detectors adjacent to one another. If two array detectors 114 
30 and 115 are used, the values are subtracted from one another to achieve dual balance 
detection. If more than two array detectors are used the signals can be selectively 
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subtracted and complex spectral density can be obtained. 

The spectral intensity as a fimction of wavelength is preferably constant. However, if 
it is not, the spectrum can be shaped in the reference, sample and/or source arms to 
make it constant Spectral shapers are known in the art. 

5 Processing 

The signal of each detector 1 14 is band pass filtered around the signal frequency, such 
as by FFT's. The signal of all detectors 114 can be combined as explained 
hereinabove to obtain the complex cross spectral density in the frequency domain. By 
Fourier transform, the complex cross spectral density can be converted to a depth 
10 profile in the tissue. Several methods to process the complex spectral density to obtain 
depth profile information are known to those skilled in the art, such as, but not limited 
to, by obtaining at least two signals with a pi/2 phase shift in the reference arm and 
then reconstructing the complex spectral density by some linear combination of the 
two signals. 

15 Following detection analog processing includes a trans impedance ampUfier, band 
pass filter, and digitization of the signal. This signal may then be converted to 
reflectivity as a fimction of depth by the Fouri^ transform opemtion. Digital 
processing includes digitization, digital band pass filtering in either the frequency 
domain or time domain (FIR or IIR fiher) and inverse Fourier transforaiation to 

20 recover the tissue reflectivity as a fimction of depth. 

System Integration 

Processing of the multiple signals may be performed using an imaging or diagnostic 
console which performs basic operations including, mathematical image 

25 reconstruction, display, data storage. Alternatively, another embodiment, shown in 
Fig. 6, shows a standalone detection and processing system 300 that may be 
connected to OCT and/or LCI systems ahready in use. In this case, the detector 302and 
digitization may be performed in the standalone unit. The input to the standalone vinit 
would be the light combined from both reference and sample arms, as previously 

30 described. The output of the system would be an interferometric signal similar to 
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previous OCT or LCI console inputs, but witii increased SNR. The standalone unit 
would contain a splitter 304 for splitting the wavelengths into spectral bands, multiple 
detectors 302, analog electronics, including TIA's 306 and an arrangement for 
reconstructing the interferometric signal, as previously described. The arrangement for 
5 reconstructing the interferometric signal would include either analog or digital 
arrangement where the analog arrangement includes band pass filters C*BPF's") 308, 
and analog arrangement for adding the individual interferograms from each 
wavelength band. The digital arrangement would include an analog to digital 
converter, and a CPU 310 capable of recombining the interferograms from each 
10 spectral band into a single friU bandwidth interferometric signal. The reconstructed 
interferogram may be then the output of ttie standalone systCTi or alternatively, the 
reconstructed interferograms demodulated signal may be used as the input to the pre- 
existing system console. 

Scan Range of the Reference arm. 

15 The ranging depth in the sample 130 is determined by the resolution with which the 
cross spectral density can be determined. In a method using a single detector the 
spectral resolution of the complex spectral density is determined by the scan range of 
the reference arm. The larger the scan range, the higher the spectral resolution and the 
larger the ranging depth in the sample. In a system with a spectral separating unit and 

20 multiple detectors, the resolution of the cross spectral density is a combination of 
reference arm scan range and spectral separating characteristics. 

Any suitable wavelength band shape may be used for separating. For arbitrary spectral 
band shapes, the scan range of the reference arm 18 is determined by the delay that is 
needed to completely resolve the spectral components in each band. 

25 For instance, in one preferred embodiment, as depicted in Fig. 7, a spectral separating 
imit can split the spectrum into two bands where each band consists of a set of narrow 
spectra in a comb-like structure. Fig. 7A shows the spectral band at detector #1. Fig. 
7B shows the spectral band at detector #2. Fig. 7C shows the combined spectral band 
of both detectors. Interleaving the comb-like spectral bands of each detector 24 gives 

30 back a continuous spectrum. The resolution needed to resolve the spectrum at an 
individual detector is half of what it would need to be in a single detector system, and 
thus the scan range of the reference arm can be reduced by a factor of two, while 
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• maintaining the same ranging depth in the sample 130. In an altemative embodiment, 
the spectral separating unit can be in the reference arm. In Fig. 8 an example is shown 
for splitting up the spectrum in several spectral bands. In this example the scan range 
of the reference arm can be reduced by a factor relating to the number of spectral 
5 bands while maintaining the same ranging depth in the sample. 

Embodiments of flie wavelength separating filter 

Several techniques are known to separate or disperse the spectrum. One method uses a 
grating and a micro lens array to focus spectral components onto individual detectors. 
A second method uses prisms instead of a grating. A third me&od uses a grating and 

10 an addressable mirror array (such as, but not limited to, a "MEMS" mirror or digital 
light processing C*DLP") apparatus or the like) to direct spectral components to 
individual detectors. A fourth method uses a linear array of optical filt^ prior to the 
array of individual detectors. A fifth method uses waveguides etched into a material or 
manufactured from fiber optic components to generate a pattern with the desired filter 

15 action. As an example, in Fig. 8 an exemplary embodiment of a wave guide filt^ is 
provided that spUts the spectrum into bands. A sixth method would use arrayed 
waveguide gratings ("AWG") to create the interleaved or arbitrary spectral bands. 

Relative intensity noise 

One of the noise terms that is present at the detectors is relative intensity noise 
20 CHIN") or Bose-Einstein noise. RIN noise likely becomes dominant over shot noise 
for spectral widths less than a few nanometers. For many detector configurations, the 
spectral width at each detector may likely be smaller than a few nanometers, and the 
relative intensity noise can dominate the overall system noise. Thus, balanced 
detection, can preferably be implemented to eliminate the RIN. Several methods 
25 known in the art exist to implement balanced detection. One such method will be 
discussed below in fiirther detail. For example, but not by way of limitation, as shown 
in Fig. 9, light from the reference arm 400 and sample arm 402 is incident on a grating 
404 at slightly dififerent angles and reflected and focused onto a linear N x M photo 
detector array 406. Along the N direction (column) of the array, wavelength is 
30 encoded. Along the M direction (row) of the array, the interference pattern of the 
sample and reference arm at a particular wavelength is recorded. Since sample and 
reference arm light were incident at slightly different angles, a pattern of interference 
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maxima and minima will be present in the column direction. Balanced detection can 
be implemented by subtracting diode signals that are exactly out of phase with respect 
to the maxima and minima pattern. Alt^atively, balanced detection can be 
implemented by measuring the amplitude of the interference pattern in the column 
5 direction which may be accomplished by subtracting the maxima or the interference 
pattem from the minima of the interference pattem along the column. An alternative 
embodiment for balanced detection is combining the reference and sample arm light 
400, 402 to produce two outputs that have interference signals with a n phase shift 
between them. This may be accompUshed by taking both output ports of a beam 

10 splitter or other beam-recombining element. The two signals may then be detected 
separately and subtracted. Since the signals that contain the interference terms are 
shifted by K in phase, these temis add constmctively upon the operation of 
subtraction. The portion of signal that contains RIN, however, cancels upon 
subtraction* The subtraction operation can occur for all M elements and be conducted 

IS in the analog or digital domain. If subtraction is performed in the analog domain, the 
bandwiddi of the signal is reduced by a factor of 2, preferably decreasing specified 
parameters of the digitization and data transfer across the computer bus. 

An example of such balanced detection is shown in Fig 10, which is described more 
fully hereinbelow. The balance detection outputs are subtracted to generate a balanced 
20 signal that cancels RIN. 

Signal processing to reconstruct the signal after spectral separatiag and detection. 

Two cases will be discussed below as nonlimiting illustrations of exemplary 
embodiments of the present invention, firstly the case of continuous spectral bands 
(blocks), and secondly the comb-like spectral bands as depicted in Fig. 7 

25 Case A: Continuous spectral bands. 

The detection arm light is split into N spectral blocks, where each spectral block 
contains the intensity between two optical frequencies, 

(20) 

The signal for the full spectral width is obtained by an FFT of the signal in each band, 
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an optional compensation of dispersion and other corrections to fixe phase and 
amplitude of each Fourier component to optimize the signal and to correct the spectral 
dmsity for side lobe reduction, addition of the complex FFT spectra, and inverse FFT 
on the added complex FFT spectrum, optionally with data reduction before the inverse 
5 FFT, to obtain the optionally demodulated function R(t), which is the interferometric 
response for a depth scan with the full source spectrum. 

Case Bl: Comb like spectral bands and the reconstruction of the full depth range in 
the sample arm from reduced reference arm scans. 

The following description provided below describes the principle of reconstruction of 
10 the full depth range in the sample arm from reduced reference arm scans according to 
the present invention. The procedure shall be explained in the case of separating the 
spectrum in two spectral bands. The exemplary method can be expanded for 
separating into many spectral bands. 

The signal at the detector for a single detector system is defined by R(t). The depth 
15 range in the sample is given by the measurement time T of a single A^line (depth 
profile) times the group velocity generated by the reference arm delay line, 
z =v T 

^ range ^ g-^ 

The smallest resolvable frequency after tan FFT is given by 1/T, which gives a 
smallest resolvable angular fi^uency Aa> = 2;r/r. The filter as depicted in Fig. 8 
20 splits the signal into two bands with peaks at m — CDQ^o)^ -hlAio^a?^ +4Afi?,efc. and 
(2? = fi>Q + Aty , o)q + 3Afi>, etc. , respectively. 

By(t) and B2(t) are the signals in band one and two respectively. The signal in 
spectral bands one and two after Fourier transform are given by 
B^ia)) = Ria))cos\cDT/4) and B2ia)) = Ria))sixx^(a)T/4). 

25 This product in the Fourier domain can also be written as a convolution in the time 
domain. Assuming the signals periodic with time T, the signals B^(t) and ^jCO are 
given by B,it)^R(t) + Rit+T/2) and B^(t)^R{t)-R(t'hT/2) . 

Using the above equations, the signal R(t) from t =0 to t=T can be reconstructed from 
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the signals (t) and B2 (/) recorded from t = 0 to t=T/2 by writing, 

i?(0 = 5,(0 + ^2(0 and i?(/ + r/2) = 51(0-^2(0 for 0<t<r/2. For higher N>2, the 
identical procedure is performed such that R(t) is reconstructed from Bj to JB^. 

This demonstrates fliat the signals B^ (t) and B2 (t) only need to be recorded over half 

5 the depth range z,^,,^^ . Thus, the depth ranging in the reference arm can be reduced by 

a factor of 2, while the ranging depth in the sample remains the same. If the signal is 
split into more spectral bands, like shown in Fig. 7, a similar procedure as described 
above allows reduction of the depfli scan in the reference arm by a factor of N, while 
the ranging depth in the sample remains the same, and N the number of spectral 
10 bands. 

An exemplary flow diagram of the procedure described above is shown in Fig. 1 1 . 
Case B2. Lumit of large number of spectral bands 

In the limit of a large numb^ of spectral bands, iV^ , the optical path length change 

in the reference arm approaches that of a wavelength, A. In this limit, only a phase 
IS change across one wavelength is needed for reconstructing flie entire axial scan over 
length L. In this case, the reference arm path delay may be accomplished by using any 
of tiie aforementioned ways for scanning the reference arm delay. Other preferred 
methods according to the present invention include insertion of an electrooptic 
modulator, acoustooptic modulator or phase control rapidly scaiming optical delay 
20 line C^RSOD") in the reference arm path to impart the path length delay of one 
wavelength. Also in this case, the wavelength separating unit does not separate the 
wavelengths into a comb pattern, but separates the spectrum into unique optical 
frequencies, with each frequency detected by a single detector. 

Case C. Fourier domain reconstruction for arbitrary wavelength patterns 

25 In contrast to the reconstruction of the LCI or OCT signal in the time or space 
domains, the signal may be reconstructed in the Fourier domain by adding the 
complex spectral components for each wavelength band to compose the Fourier 
transform of the LCI or OCT signal. Alterations of the phase for each Fourier 
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component may be preferred in certain selected circumstances to correct for 
nunimization of reference aim delay length. 

Reconstruction of the image or one dimensional axial scan 

5 Following reconstruction of the LCI or OCT signal in the real domain, the axial 
reflectivity may be determined by demodulating the reconstructed LCI or OCT signal. 
An arrangement for demodulation can include multiplication by a sinusoid and low 
pass filtering, envelope demodulation using envelope detection, square law 
demodulation and low pass filtering, quadrature demodulation followed by FIR, IIR 

10 filtering, or low pass filtering. In addition, the reconstruction of Stokes vectors 
(polarization) and flow from these LCI or OCT signals is known to those skilled in the 
art. Following reconstruction and demodulation, the data may be displayed in one or 
two-dimensional format (image) for interpretation and ultimately diagnosis of a tissue 
condition or defect in a medium. If the LCI or OCT signal is reconstructed in the 

15 Fourier domain, such reconstructed signal in the Fourier domain can be demodulated 
in the Fourier domain by shifting tiie Fourier spectrum and pfenning an inverse 
Fourier transform. As a result, the complex signal in the real domain (quadrature 
signal) is then reconstmcted into axial reflectivity information by computing the 
amplitude of the real portion of the quadrature signal. The complex component is used 

20 for computing polarization or flow information. Alternatively, if the signal is 
reconstructed in the Fourier domain, it can be directly inverse Fourier transformed 
into the real domain and undergo the aforementioned processing described for the 
reconstructed real domain signals. 

Fig. 12 shows an exemplary embodiment of a Spectral Domain OCT interferometer 
25 design 500 showing spectral compounding of light sources 502, 504, and 506 and 
acousto-optic generation of the carrier in tiie reference arm. The blocks labeled AOM 
are acousto-optic modulators 508, 510. The two outputs each go to separate spectral 
detection units 1 14, 1 15 (as depicted in Figs. 3 and 13) for balanced detection. 

After spectral compounding of the source li^t in the first 50/50 splitter and the 80/20 
30 splitter, light enters a modified Michelson interferometer. A configuration that 
implements balanced detection is shown. The sample arm goes to the probe (e.g., a slit 
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lamp). Reference ann light is transmitted flarough two acousto-optic modulators with 
a difference frequency of 10 kHz to generate a constant carrier frequency that is 
independent of wavelength. The balanced detection outputs go to separate spectral 
detection imits. 

5 Spectral detection unit 

Referring to Fig. 13, the core of Spectral Domain OCT is spectral separation of the 
detection aim ligjit onto a multi-element array 114. The detection arm beam 520 is 
spectrally separated by a grating 520 and focused by a lens 522 onto a multi-element 
array 114. 

10 A scan cameras with N detector elements is used as spectral detection unit 128 (see 
Fig. 2). Preferably, balanced detection is implemented by adding a second line scan 
camera. As is known to those skilled in the art, the depth range is inversely 
proportional to the spectral resolution. When the real part of the complex spectral 
density is determined, ranging depth z is defined by, 

15 ^=4h' (i^> 

4A1 

Line scan rates of 20 kHz can be achieved, allowing demodulation of a 10 kHz carrier 
to extract the complex cross-spectral density. Data is digitized and transfCTred to 
computer memory. Demodulation of the signal is done in software. Scan rates of 
10,000 depth profiles per second or more can be achieved. 

20 Dual Balanced detection 

Dual balanced detection is preferably used by the present invention, which is 
preferably utilized for the following reasons. Firstly, most ligjit sources generate 1/f 
noise (f = frequency) at relatively low frequencies. Balanced detection will eliminate 
1/f source noise. Secondly, an interference term of the sample arm light with itself 
25 (auto-correlation term) is present on top of the true signal term, which is the 
interference between sample and reference arm. This auto-correlation term can be 
eliminated by a differential technique. Balanced detection may eliminate this auto- 
correlation term from the measured signal. Thirdly, RIN can be reduced. 
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Data acqiiisition and processing unit 

The data rate at 20,000 spectral profiles per second, with 2000 detector elements and 
8-10 bit resolution (the dynamic range of most line scan cameras) is 40-80 MB/sec. 
Maximum sustainable data transfer speed over the PCI bus is 100 MB/sec. In a 
5 computer with two independent PCI bridges to computer system memory, 
approximately 200MB/sec of data can be transferred for real time processing of data 
fi'om two line scan cameras simultaneously. Implementation of dual balanced 
detection in analog by subtracting line scan camera signals before digitization may 
reduce the data rate by a factor of 2. High-speed data acquisition boards are available 

10 at resolutions of 12-14 bits and speeds up to 100 Msamples/sec. A single 2048 point 
fast Fourier transform on a 2.5 GHz Pentixmi 4 processor takes 50 jisec. These 
numbers show that real-time processing of Spectral Domain OCT data at 20,000 
spectral profiles/sec is within reach of current data acquisition and processing power 
of dual processor PC's. The data collected by the spectrometer can be sampled with 

15 equal wavelength increments. Fourier transform, however, links z and k space (or t 
and w). Because of the non-linear relation between k and X the spectrum jfirom the 
spectrometer should be interpolated to create evenly spaced samples in k domain. To 
achieve the optimal point spread function, dispersion in the sample and reference arm 
of the interferometer should be balanced. We have shown that dispersion imbalance 

20 can be corrected for by digital processing, allowing for correct compensation of 
dispersion for individual eye lengths. 

PHASE TRACKING 

The present invention also provides apparatus and methods for phase tracking in 
spectral domain ("SD") OCT. 

25 FuUy paraUel SD OCT 

One of the features of fully parallel SD OCT is spectral dispersion of the detection 
arm light onto a multi-element array such as but not limiting to an integrating device 
(e.g., CCD) and measurement of the real or complex spectral density at high speeds. 
The detection arm beam is separated by a spectral separating unit (e.g,, grating) and 
30 focused onto the array. With respect to previous Spectral Domain OCT designs 
known in the art, two differences are apparent that will be discussed below: 1) 
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implemratation of balanced detection, and, 2) implementation of phase tracking. 

Spectrometer design The depth range in SD OCT is inversely proportional to the 
spectral resolution. Using the complex spectral density, ranging depth z is given by, 

2/iAA (Ig^ 

5 Dual balanced detection: Dual balanced detection is advantageous for at least three 
reasons. First, most light sources generate 1/f noise at relatively low frequencies (tens 
of kHz range). In time domain CTD") OCT systems 1/f noise is not a problem 
because the signal carrier is in general in the MHz range where 1/f noise is not 
significant In SD OCT, balanced detection may likely eliminate 1/f source noise. 
10 Second, an interference of the sample arm light with itself (auto-correlation term) is 
preset on top of the true signal. This auto-correlation term can be eliminated by a 
differential technique. Balanced detection can be used to eliminate this auto- 
correlation lean from the measured signal. Third, balanced detection may reduce 
relative intensity or Bose Einstein noise. 

IS Phase Tracking: Phase tracking is preferable to eliminate phase instabilities in the 
interferometer. Phase instabilities can cause individual interferometric fringes to shift 
in location. If detection is slow relative to the shifting of the fringes, the resulting 
averaging results in an artifactual decrease in the measured fringe amplitude. Fast 
detection arrays can capture the cross spectral density at a rate of 20 to 40 kHz, 

20 resulting in integration times of 50 to 25 psec, respectively. Phase instabilities arising 
on a time frame shorter than the integration time of tiie array should be compensated. 

Fig. 14 shows an exemplary interference pattem as a function of path length 
difference between sample and reference arm. 

Phase locking circuitry is common in electronics, and is frequently used in radar and 
25 ultrasound. Active phase tracking can be implemented by modulating the 
interferometer path length difference at 10 MHz with an electro-optic phase modulator 
in the refi^ence arm over a fraction of the wavelength. By demodulating the intensity 
measured by one detector at the output of the interferometer at the frequency of the 
path length modulation, an error signal can be generated indicating in which direction 
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the phase modulator should shift to lock onto a fringe amplitude maximum. By adding 
an offset to the phase modulator as determined by the error signal, the phase tracker 
actively locks onto a fringe maximum. The phase modulator can only modulate the 
path length difference over a few wavelengths. The processing unit can determine if 
5 the phase modulator has reached its range limit, and jump by a full wave in phase to 
maintain lock on a different fringe maximum. This approach exploits the fact that 
phase should be controlled only modulo 2n. In addition, the processing drives a 
slower component (e.g., the Rapid Scanning Optical Delay line) to extend the path 
length range of the phase modulator/RSOD combination over several millimeters. 
10 Phase locking can be performed on a fringe maximum, minimum, or zero crossing, 
based on the type of mixing performed in the demodulation circuit. 

The present invention can also use autoranging technology, including processing 
algorithms, as disclosed in copending U.S. application No. 10/136,813, filed April 30, 
2002, entitled METHOD AND APPARATUS FOR IMPROVING IMAGE 
15 CLARITY AND SENSITIVITY IN OPTICAL COHERENCE TOMOGRAPHY 
USING DYNAMIC FEEDBACK TO CONTROL FOCAL PROPERTIES AND 
COHERENCE GATING, and commonly assigned to the assignee of the present 
invention, the disclosure of which is incorporated herein. 

The autoranging mechanism may, in one exemplary embodiment, comprise a 
20 processor imit for (a) obtaining a first scan Une; (b) locating a surface location "S" of 
a sample; (c) locating an optimal scan range '"R" of the sample;(d) modifying a 
reference arm delay waveform to provide an output; (e) outputting the output to a 
reference arm; (f) determining whether the image is complete; and (g) moving to the 
next scan line if the image is not complete or remapping the image using the surface S 
25 data and the waveform data stored in the memory storage device if the image is 
complete. 

If the light retumed from the sample is of low ampUtude, phase locking may be 
xmstable due to the presence of noise. In another embodiment, a separate, preferably 
monochromatic, light source is input into the interferometer. The separate source 
30 wavelength may overlap with the broad bandwidth OCT or LCI source spectrum or 
may be centered at a different wavelength than the OCT or LCI source spectrum. The 
separate source is preferably of higher power and may be combined with the source 
arm (using wavelength division multiplexer, grating, prism, filter or the like) travel to 
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the reference and sample arms and return back to the beam recombining element. The 
returned separate source light can then separated from the OCT or LCI light following 
transmission back through the beam recombining element (i.e. beam splitter output). 
A separation arrangement can perform spectral separation by a dispersing element, 

5 such as a dichroic mirror, filter, grating, prism, wavelength division multiplexer or the 
like. The separate source will be detected separately from the OCT or LCI broad 
bandwidth light using one or more detectors. The higher power provided by this 
separate source can enable detection of a higher amplitude interference pattern, and 
provide an improved input to the phase tracker, thus enabling more stable phase 

10 tracking. 

Fig. 15 shows one exemplary embodiment of a phase tracker system 600 according to 
the present invention with an extended phase lock range, by combining a fest element 
(EO phase modulator) 602 to modulate the path length difference over a small range, 
and a slower element (RSOD) 604 to modulate the path length over an extended 

IS range. The detector 606 signal is mixed with the phase modulator modulation 
frequency 608 by a mixer 610 and low pass filtered (filter not shown) to generate an 
error signal. The processing unit 612 preferably processes the error signal to generate 
an offset voltage, and adds this offset voltage to the modulation signal 608, so as to 
generate the output for the phase modulator driver 614. In addition, the processing 

20 unit 612 can generate a signal to the RSOD 604 to provide extended range tracking of 
the phase over distances of several millimeters. Light source 616, fiber splitter 618, 
sample arm 620 and reference arm 622 are shown, and are described herein. 

Mixer Implementatioii: The intensity I(t) at the detector at a given moment within a 
single oscillation of the fiinge pattern is given by 

25 /(O = cos[^(0] 

where the phase (p gives the position in the fringe. For (p=0, the signal is at a fringe 
maximum, for (p=7c, the signal is at a fringe minimum. At an arbitrary moment t, the 
phase 9(t) is given by, 

30 where a desoibds the position within a single oscillation of the fringe pattern, and 



28 



wo 03/062802 



PCT/US03/02349 



p*sin(a)t) is the phase modulation introduced by the phase modulator, with p the 
amplitude of the phase modulation, and o) the frequency of the phase modulation 
signal. The intensity at the photodetector I(t) can be mixed with a carrier at frequency 
(0 and 2co, resulting in the mixer signal MixerC(t), MixerS(t), Mixer2a)C(t) and 
5 Mixer2o)S(t), 

MixerC(t) = cos(^/) * cos(a + sin(^O) • MixerS(t) = sin(6;/)* cos(a + sm(a)t)) 

MixerlQ)C{t) = cos(2cyr) * cos(a + J3 sin(fi;0)- 
Mixerla^it) = sin(2<y^)* cos(ar + /3 smicot)) 

The time average over a single oscillation of the carrier frequency co of MixerC, 
10 MixerS, Mixer2coC and Mixer2a)S is given by, 

MixerC(t) = 0 . MixerSit) = sin(a) * J^{/3) . MixerlwCit) = cos(a) * J^^ifi) . 

3 5 » 

Mixer2(oS(f) = 0 

where J,(P) and JsCp) are a Bessel frmctions of the first kind; its value depends on p, 

the amplitude of the phase modulation. Thus, the signal ^^^^(0 and J^ixerlcoCit) 
15 are proportional to sin(a) and cos(a), respectively, with a the position within a single 

oscillation of the fringe pattern. The mixer outputs ^^^^^it) Mixer2(oC{t) 
used as an error signal to generate an offset voltage to steer the phase modulator to a 
new center position that minimizes the error signal, and locks the interferometer 
output on a fringe maximum or minimum, or a zero crossing, respectively. The 

20 complex spectral density can now be determined by two consecutive array scans, one 
where the error signal sin(a) is minimized, and the next where the error signal cos(a) 
is minimized, resulting in a 90 degrees phase shift between the two interference 
patterns. Using this mixing arrangement, the complex spectral density can be obtained 
rapidly and without resorting to an additional mechanical arrangement for changing 

25 the phase of the reference arm light. 

Fig. 10 shows one exanplary embodiment of a SD OCT system 700 with phase 
tracker for providing balanced detection according to the present invention. In this 
embodiment, a source 702 provides light which passes through a splitter 704, which 
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sends part of fhe light to a sample probe 706 and the remainder of the light to a Rapid 
Scanning Optical Delay C^RSOD'O line 708. Light is passed fix)m the RSOD 708 to 
the phase modulator PM 710. Light from the phase modulator PM 710 is sent through 
a splitter 712, and then through two additional splitters 714 and 716, a portion of the 
5 output of which is sent as balanced detection outputs to spectral detection units (not 
shown, but as described elsewhere herein) and the remainder of the output is sent to 
the phase tracker assembly 720. In the phase tracker assembly 720, phase tracker 
detectors and Dj, 722 and 724, receive the partial output of the pair of splitters 714 
and 716, which in turn send signal to a mixer 726 to generate an error signal. A 

10 processing unit 728 processes the error signal, where the sum generation of offset 
voltage and adds this to the modulation signal 730 to generate the output for the phase 
modulator driver 732. Modulation signal, shown at box 730, is forwarded to the mixer 
726 and the processing unit 726. In addition, the fringe amplitude could be too small 
for the phase tracks: to lock. Alternatively, a secondary source with longer coherence 

IS length could be coupled to the system 700 to provide a larger fringe amplitude to the 
phase tracker. 

The present invention provides a method for tracking phase in an imaging system, as 
shown in Figs. ISA-C the method comprising the steps of: (a) measuring a signal 
received from the sample arm; (b) increasing a phase of the signal; (c) measuring a 
20 first signal partition of the signal defined as x^ at at least one peak of the signal; (d) 
deteraiining whether to increase or decease the phase of the signal by an incremental 
amount; (e) after step (d), measuring a second signal partition of the signal following 
step d); and, if the signal is at its peak, remeasuring the signal and if the signal is not 
at its peak, repeating steps d) and e). 

25 The method fiarfher may comprise that steps (a) - (f) are performed in parallel with 
o&er imaging processes. The adjustment of phase "9" is defined as A(x2-x,), where 
**A" is a constant. Furthermore, optionally, step d) may frirther comprise the substeps 
of dl) determining wheth^ A(x2-x,) is within range of the phase modulator; and d2) 
changing <p by an amount equal to ACxj-x^) if A(x2-x,) is within the range or changing 

30 q> by an amount equal to A(x2-Xi)-m27i: if ACxj-Xj) is outside of the range, where M is 
an integer greater than 1. The method may optionally ftirflier comprise a substep d3) 
remeasuring signal x,. 
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Data acquisition and processing unit 

In general, the data collected by the spectrometer are sampled with equal wavelength 
increments. Fourier transfonn, however, links z and k space (or t and w). Because of 
the non-linear relation between k and X the acquired spectrum is interpolated to create 
5 evenly spaced samples in the k domain. Alternatively, the light could be dispersed in 
such a way on the detection array that the light is samples in equal intervals in k 
space, such that the interpolation becomes obsolete. Alternatively, the detection array 
spacing could be designed to sample the light evenly spread in the k domain, such that 
the interpolation becomes obsolete. To achieve the optimal point spread function, 
10 dispersion in the sample and reference arm of the interferometer should preferably be 
balanced. Dispersion imbalance can be corrected by digital processing. 

The present invention provides a probe for locating atherosclerotic plaque in a blood 
vessel, comprising: an interferometer; a spectral separating unit which splits signal 
received from the interferometer into a plurality of optical frequencies; and a detector 
15 arrangement capable of detecting at least a portion of the optical frequencies received 
from the spectral separating miit. 

The present invention fruHier provides an apparatus for delivering a therapeutic agent, 
comprising: a probe disposed in the housing and comprising: an interferometer, a 
spectral separating unit which splits signal received from the interferometer into a 

20 plurality of optical frequendes, a detector arrangement capable of detecting at least a 
portion of the optical frequencies received from the ^ectral separating unit; and a 
conduit cooperating with the probe, and comprising a proximal end for receiving the 
therapeutic agent and a distal end for delivering the therapeutic agent at a 
predetermined location, the location being determined by imaging the environment in 

25 proximity to the distal end using the probe. 

An exemplary embodiment of the present invention will be ftirther described below in 
connection with the following example, which is set forth for purposes of illustration 
only. 

EXAMPLE 

30 The method according to the present invention was verified in the laboratory by the 
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following e?q)eriment. 

In the existing OCT system, the shot noise power spectrum as detemiined from the 
spectral density due to the reference arm optical power was measured. Then 2/3 of the 
spectrum from the referCTice arm was blocked, and experimentally it was verified that 
5 the shot noise power spectrum was reduced by a factor of three, thus demonstrating 
that the shot noise is reduced by a factor of 3 if the spectrum is split in three spectral 
bands (see Fig. 16). The upper curve (gray dotted line) shows the power spectrum for 
the OCT signal with one detector. For the lower curve (solid line), the spectrum was 
limited by 1/3 with a corresponding factor of 3 improvement in signal to noise ratio. 
10 This data was generated by experiment, blocking 2/3 of the spectrum in a grating- 
based double-passed pulse shaping rapidly scanning optical delay line. 

An object with low reflectivityr was inserted in the sample aim. Using the ftiU spectral 
width of the source, the power ^ctrum of the interference between sample and 
reference arm light was determined in the lower half of the spectral density. Then the 

15 upper part of ttie source spectrum was blocked in the reference arm, and it was 
verified that the lower 1/3 of the power spectrum of the interference between sample 
and reference arm light had the same magnitude as in the previous measurement (see 
Fig. 17). This figure demonstrates that the signal amplitude is equal for the N=l and 
N=l/3 cases where they overlap. The result of equal amplitude signal for N=l/3 case 

20 and the 3-fold lower noise for the N=l/3 case (see Fig. 6) demonstrates that splitting 
into N wavelength bands increases the SNR by a falctor of N. 

This demonstrates that when the light in flie detection arm is split in two spectral 
bands, the spectral density of flie interference between sample and reference arm light 
within the spectral bandwidth of a single detector is imchanged. Combined with the 
25 measurement that showed a reduction in the shot noise power spectrum, the 
conclusion is that a reduction of shot noise can be realized by splitting the detection 
arm light in separate spectral bands. 

Experimental verification of the noise reduction. 

To demonstrate the noise reduction in Spectral Domain OCT, an OCT system was 
30 used, including a Rapid Scanning Optical Delay line (RSOD) was used in the 
reference arm, enablmg portions of the spectrum to be blocked. Detector signals were 



32 



wo 03/062802 PCTAJS03/02349 



digitized at 2.5 Msamples/sec, allowing digital processing of the fi±tige information. 
First, the thermal noise density of the detector was measured as a function of 
frequency by blocking all light onto the detector. Second, the shot noise density of the 
reference arm power was measured with only the reference arm power incident on the 
5 detector. Third, both the sample and reference arm light were incident on the detector. 
The sample was a single scattering surface mounted in a model eye and 512 depth 
profiles were acquired in 2 seconds. The power density I{ff was measured, which 
is proportional to the spectral density squared (see Eq. (9)). Then we blocked half of 
the spectrum in the reference and measured again the shot noise density of the 
10 reference arm by blocking the sample arm, and the power density I{fY when both 
sample and reference arm light were incident on the detector. Shot noise and power 
densities were corrected for thermal noise by subtraction. Thermal noise was at least a 
factor of 3 smaller than the lowest shot noise level. 

Fig. 18 shows a graph of power density for the full spectrum, and with half of the 
15 spectrum blocked in the reference arm, as a function of frequency. The solid line 
shows the power density for the full spectrum. The shot noise level measured while 
the sample arm was blocked is also shown. The dashed line shows the power density 
with half the spectrum blocked in the sample arm. The shot noise level measured 
while the sample arm was blocked is also shown. Fig. 18 demonstrates that the shot 
20 noise level was reduced by a factor of 2 by blocking half the spectrum in the reference 
arm. At the same time, the signal at frequencies corresponding to wavelengths that 
were not blocked in the reference arm remained the same. 

As is evident from Fig. 18, which sununarizes the measured results, the shot noise 
daisity is reduced by approximately a factor of 2 by blocking half the spectrum in the 

25 reference arm. Fig. 19 shows that after subtraction of the shot noise levels from the 
corresponding signals, the power densities for those frequencies that corresponded to 
wavelengths that were not blocked in the reference arm remained the same. This 
demonstrates that the shot noise density is reduced by a factor of 2 when the total 
reference arm power is reduced by a factor of 2 by blocking half the spectrum, while 

30 the signal power density for wave lengths not blocked in the reference arm remains 
unchanged. 

Fig. 19 shows a graph of the square root of the power densities for the ftiU spectrum, 
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and for half the spectrum blocked in the reference aim as a function of frequency. The 
solid line shows the spectrum after subtraction of the respective shot noise. The 
dashed line shows the half spectrum after subtraction of the respective shot noise. Fig. 
13 demonstrates that after subtracting the respective shot noise contributions, the 
5 signal at frequencies corresponding to wave lengths that were not blocked in the 
reference arm remained the same. 

The next experiment fiirther demonstrated that by dispersing the spectrum in the 
detection arm over several detectors, and by selectively band pass filtering the signals 
of each detector, the SNR is increased. The detection arm light was dispersed over 4 
10 detectors by a diffraction grating as shown in Fig. 13, and the detector signals were 
separately amplified by transimpedance amplifiers with a bandwidth of 600 kHz and 
simultaneously digitized. 

Fig. 13 show;5 a schematic of an exemplary apparatus setup used to demonstrate SNR 
improvement by Spectral Domain OCT according to the present invention. Scanning 
15 of the reference arm 106 was performed witib a Rapid Scanning Optical Delay line 
(RSOD) 120. Individual signals from the array detector 114 were amplified by 
transimpedence amplifiers, digitized by a 4-charmel 2.5 MHz per channel A/D board 
and stored in computer memory (not shown). 

First, the th^mal noise density of all four detectors was measured. Second, the shot 
20 noise density of the reference arm light in each detector channel 1 16 was measured. 
Third, both the sample and reference arm light were incident on the detector 1 14. The 
sample 130 was a single scattering surface mounted in a model eye and 512 depth 
profiles were acquired in 2 seconds. The power density /(/)^ in each detector chaimel 
1 14 was measured. Then, the signals of the four detectors 114 were summed, and the 
25 combined power density /(/)^ was determined. The results are shown in Fig. 20, 
which demonstrates fliat the shot noise is lower in each individual channel compared 
with the sum of all channels, but that the power densities /(/)^ in the individual 
chaxmels within their respective bandwidths are approximately equal to the power 
density I{fY of the coherent sum of the four chaimels. 

30 Fig. 20 shows a graph of the power densities for four separate detectors 1 16 of Fig. 
13. The spectrum in the detection arm was dispersed over four separate detectors 116 
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by a difiBraction grating 520. The shot noise levels for each individual detector 1 16 are 
significantly lower than fot the coherent sum of the four detector channels. Bars at the 
top of the image indicate the signal pass band that was applied to the individual 
channels and the coherently summed channel to generate Fig. 21. 

5 In Fig. 21, the square of ttie coherence envelope is shown for both the direct sum of 
all four detection channels and the coherent sum after digitally band pass filtering 
each detector channel with a bandwidth centered at the center firequency of the 
respective detector signal. Fig. 21 shows that the interference fringe signal I(t) of the 
direct sum and the band pass filtered coherent sum of the four detector signals results 

10 in virtually the same coherence envelope peak value, while the band pass filtered 
coherent sum of the four detector signals shows a significantly lower noise level than 
the direct coherent sum. Since the pass band of each individual channel was slightly 
larger than one third of the pass band of the fiill signal (pass bands are indicated in 
Fig. 20), an increase of SNR of a factor of 2.87 was expected. The noise level dropped 

15 by a factor of 2.8. However, band pass filtering also reduced tihie signals slightly, by a 
factor of 1 , 12, resulting in an effective increase in SNR of a factor of 2.5. 

These experiments clearly demonstrate that spectrally dispersing the light in the 
detection ami can offer a significant SNR advantage. 

Fig. 21 shows a plot of the coherence envelope for the coherently summed chamels, 
20 and the coherently sunamed channels after band pass filtering each channel. The solid 
line is the sum of channels. The dashed line is the pass filtered sum of chaimels. Fig, 
21 clearly demonstrates the signal to noise gain that can be achieved by spectrally 
dispersing the signal in the detection arm over several individual detectors. In this 
example the noise level was reduced by a factor of approximately 2.8. Since the 
25 coherence peak was reduced by a factor of 1 . 12 due to some remaining signal firaction 
filtered out by the band pass filters, the actual SNR improvement was 2.5. 

Although only a few exemplary embodiments of this invention have been described in 
detail above, those skilled in the art will readily appreciate that many modifications 
are possible in the exemplary embodiments without materially departing from the 
30 novel teachings and advantages of this invention. Accordingly, all such modifications 
are intended to be included within the scope of this invention as defined in the 
following claims. It should further be noted that any patents, applications and 
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